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1 Introduction 

In the last two chapters, we would like to model asset dynamics and price options (or other deriva- 
tives). Recall that we are able to price a forward contract without assuming any model. This is 
not the case for option valuation. In this course, we will discuss two ways to model asset dynamics: 

(1) Binomial Lattice Model (Discrete Time) 

(2) Brownian Motion and Stochastic Calculus (Continuous Time) 


2 Binomial Lattice Model 


Setting: 

• S = Sq is known. 

• After one period (At), S becomes uS (upward) or dS (downward) 

• u > 1 and 0 < d < 1. 

• Probability of going up is p. 

• Assume the following two quantities are known: 

(1) Expected annual rate of return: v = E[log(St+i/<St)]. 

(2) Variance of the annual rate of return: a 2 = Varflog^+i/St)]- 

Using 


We find that 


vAt. = E(\og(S A t)) 
cr 2 At = Var(log(SAi)) 

1 


1 1 1 


u = ex pji/^At) 2 + cr 2 Af} 
d = exp{ — \J (vA t) 2 + a 2 At}. 
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3 LOGNORMAL DISTRIBUTION 


For At small, we can use 


P = 


u 

d 


1 1 V r~z — 

— H V At, 

2 2 a 


_ e <r\/Af 

_ g — <7%At 


Example 2.1. (2012-2013 Final Q4) Consider the price of a stock at time t = 0 given by Sq = 16.0, 
with the mean and standard deviation of the annual rate of return given by is = 10% and a = 25% 
respectively. Suppose that the risk-free rate is given by 5%. 

(a) Based on At = 1/12 (i.e. one month), draw a 3— step binomial lattice for the next three months 
Si, S 2 and S 3 . What are the probabilities of obtaining each of the possible outcomes for S 3 ?. 

(b) Estimate the price of a European Call option expiring in 2 months with strike K = 16.0 based 
on a binomial model. 

Solution, (a) u = 1.0748; d = 0.9304; p = 0.5722. 

S 0 = 16; uS = 17.197; dS = 14.886; u 2 S = 18.484; d 2 S = 13.850; u 3 S = 19.868; d 3 S = 12.885. 
P(S 3 = u 3 S) = 0.1873; P(S 3 = uS) = 0.4202; P(S 3 = dS) = 0.3142; P(S 3 = d 3 S) = 0.0783. 

(b) p* = 0.51095. 

/ = 0.6431. 


3 Lognormal Distribution 


Definition 3.1. If X ~ N(p,(j 2 ), then the random variable Y := e x is lognormally distributed. 
It is clear that Y > 0. 

Example 3.1. (2012-2013 Final Q3 Modified) Consider a random variable Y which is log-normally 
distributed, with log Y = X N(p,a 2 ). 

(a) Show that the moment generating function Mx(t) of X is exp(^t + a 2 t 2 / 2). 

(b) Show that the variance of Y is Var(E) = exp(2/i + <7 2 )[exp((j 2 ) — 1], 

(c) Compute the skewness of Y. k 3 = E([E — /iy] 3 )/crf.-. Does it suggest that the distribution is 
right-skewed? 

Remark 3.1. A distribution is right-skewed if the right tail is longer or fatter than the left side. 

Solution, (a) To find the moment generating function or deal with normal density, we often use 
the completing square trick: 

(x - p) 2 


E(e 


Xt\ 


/ 


OO 1 

Xt 


— OO 
OO 


v / 2 tt 


exp 


1 


/ — 00 \! 2 tt 

f°° 1 


/ — 00 \/ 2 tt 


exp 

exp 

1 


2 ^ ' dx 

x 2 — 2 xp + p 2 — 2 xa 2 t 

2 ^ 


dx 


[x — (p + (J 2 t )] 2 — 2 pa 2 t — aH 2 


2 a 2 


dx 


exp <j pt + -cr 2 t 2 


exp <j pt + - a 2 t 2 


I — 00 's/ 27T 


exp 


[x — (p + a 2 t)] x 

2 ^ 


dx 
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4 STOCHASTIC PROCESSES AND MARTINGALES 


(b) 


Var(Y) = Var(e A ) 

= E(e x ) 2 - [E(e x )] 2 
= E(e 2X ) - [E(e x )] 2 

— g2/.t+2<r 2 _ |- e ^+CT 2 /2j2 

_ g2^+2cr 2 _ g 2 iM+cr 2 

= e 2 ^ +a \e a2 - 1 ). 


E (Y - 


o'- 


Y 


EY 3 - 3/iyEY 2 +34EY- n\ 
(exp(2/i + cr 2 )[exp(<r 2 ) — l]) 3 / 2 


exp(3^i + |<T 2 ) — 3 exp(3/i + § o 2 ) + 3 exp(3/u + fo -2 ) — exp(3 /u + |cj 2 ) 


ex' 


p(3/x+ ^)(exp(a 2 ) - l) 3 / 2 


> 


exp(3<r 2 ) — 3exp(a 2 ) + 2 
(exp(cr 2 ) — l) 3 / 2 

(exp(er 2 )) 3 - 3exp(a 2 ) + 2 / 2 1 

(exp (o 2 )) 2 — 2exp(cr 2 ) + 1 
(exp(cr 2 ) + 2)y / exp((j 2 ) - 1 
0 . 


Therefore, log-normal distribution is right-skewed. 

Example 3.2 (09-10 Final Q3). A simple gross return, (1 + R) is lognormal (0, (0.1) 2 ), which 
means that log(l + R) is N(0, (0.1) 2 ). Evaluate P(1 + R < 0.9). You may leave your answer in 
terms of the CDF of a standard normal random variable. 

Solution. 


P(1 + i? < 0.9) = 


where is the CDF of a standard normal random variable. 

4 Stochastic Processes and Martingales 

A stochastic process Xt (or X(uj,t)) is a collection of random variables indexed by an index 
set T. It is a function from 0 x T — > M, where D is the sample space. Usually, T = {0, 1,2,.. .} 
(which corresponds to a discrete stochastic process) or T = [0, 00 ) (which corresponds to a 

continuous stochastic process). 

We can think of t as time and each a; is a particle. With this picture in mind, Xt(uj) would 
represent the position at time t of the particle (experiment) oj. 


P(log(l + R) < log(0.9)) 
p( log(l + R) < log(0.9) , 


d>( 


0.1 
log(0.9) 
0.1 


0.1 ' 

) = $(-1.0536), 
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5 WIENER PROCESS 


Definition 4.1. A stochastic process {X n ;n 6 N} is a martingale if for each n, 

(a) E(|X n |) < oo 

(h) E(X n \X 1 ,...,X n _ 1 ) = X n _ 1 

Remark 4.1. If we think of X n — X n _i as the net winnings per unit stake in game n in a series 
of game, played at time n = 1,2,.. .. Then, we have E(X n — X n _i \ X\, . . . , X n _i) = 0. So that 
the game series is fair. 

Example 4.1. Let X n be independent random variables with E(X n ) = 0 for all n > 1, then 
S n = X ! + ••• + X n is a martingale. 

Solution. 


E (S n | Si,...,S n -i) 


E(5 n _!+X n | 5i,...,5 n _!) 
Sn- 1 + E (X n ) 

S n - 1- 


Example 4.2. Let X n be a symmetric random walk, that is X n = Z\ + • • • + Z n where Z % are iid 
random variable with P(Z n = 1) = P(Z n = —1) = 1/2. Show that X/ — n is a martingale. 

(Hint: Use the properties of conditional expectation.) 

Solution. 


E(X ^ +1 — (n + 1) | X 1; ... . ,X n ) 


E((X n + Z n+l f - (n + 1) | X u . . . , X n ) 

E(X^ + 2X n X n+ i + Z^+i — (n + 1) | X \ , . . . , X n ) 
X 2 n + 2X n E(Z n _|_i) + E (Z 2 +1 ) — (n + 1) 

X 2 + 0 + 1 — (n + 1) 

X 2 — n. 


5 Wiener Process 

Definition 5.1. A stochastic process W(t ) is called a standard Brownian motion or Wiener process 

if 

(1) With probability one, !U(0) = 0. 

(2) (W(t) has independent increments) For any t\ <t 2 <t^ < 1 4 , W(t/f) — W(t 3 ) is independent 
with W{t 2 ) - W(h). 

(3) For s < t, W(t ) - W(s) ~ N(0,t- s ). 

/// Almost all sample paths ofWt are continuous functions. 

Properties: 

• W(t) is normally distributed for any t. 

• For s < t, 

E (W t | W s ) = W s . 

• W{t ) is nowhere differentiable (with probability one). 
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6 APPENDIX 


Remark 5.1. A formal introduction of Brownian motion, conditional expectation, martingale, Ito 
integral and stochastic differential equation will be given in RMSC4005 Stochastic Calculus. 

Example 5.1. Let Wf be a Brownian motion. Show that for s < t, W ( — t satisfies 

E (Wf -t\W s ) = W 2 - s. 

What can you conclude for W 2 — t? 

Solution. 

E(W 2 -t\W s ) = E[(W t -W s + W s ) 2 \W s ]-t 

= E [(W t - W s ) 2 + 2W s (W t - W s ) + W 2 \W S ] - t 
= E[(Wt-W s j 2 ] + 2W s E{Wt- W s ) + W 2 -t 
= (t — s) + W 2 — t 
= W 2 -s. 

Therefore, W 2 — t is a martingale. 

Example 5.2. Prove that TL^ 3 — 3 tWt is a martingale. 

Solution. 

E (W? - 3 tW t | W a ) 

= E[(W t -W s + W s ) 3 -3t(W t - Ws + W a )\W a ] 

= E [(W t - W s ) 3 + 3 (Wt - ng 2 w s + 3 (W t - W S )W 2 + W 3 - 3 t(W t - W S )\W S ] - 3 tW s 
= E (W t - W s ) 3 + 3W S E(W t - W s ) 2 + 3 IT 2 E(W t - W s ) + W 3 - 3 1 E(W t - W s ) - 3 tW s 
= 3W s {t-s) + W 3 -3tW s 
= W 3 -3 sW s . 

Remark 5.2. Some properties of conditional expectations: 

(1) (Independence) E(g(X)|y) = E(g(X)) if X and Y are independent. 

(2) (Taking out what is known) E(g(X)h{Y)\Y) = h(Y)E(g(X)\Y). 

6 Appendix 

This appendix summarizes some of the facts about option you learned in RMSC2001. 

6.1 Options 

Definition 6.1. An option is a contract which gives the buyer the right (but not the obligation) 
to buy or sell an underlying asset at a specified strike price on or before a specified date. 

Two basic types: 

(1) A call option gives the holder of the option the right to buy an asset by a certain date for a 
certain price. 

(2) A put option gives the holder the right to sell an asset by a certain date for a certain price. 
Terms: 
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6.2 Risk-Neutral Valuation 


6 APPENDIX 


(1) Maturity date: the date specified in the contract 

(2) Exercise price/ Strike price: the price specified in the contract 
Payoff of the options: (Note that payoff / profit. Profit is payoff minus cost) 

(1) Call option: max(iSr — K, 0). 

(2) Put option: max (A' — St. 0). 

where T is the maturity date and K is the strike price. According to when you can exercise: 

(1) European option: an option that may only be exercised on expiration 

(2) American option: an option that may be exercised on any trading day on or before expiry. 

(3) Bermudan option: an option that may be exercised only on specified dates on or before 
expiration. (Bermuda is between Europe and America!) 

If the structure of the option is different from standard calls and puts (plain vanilla products), 
then it is referred to as an exotic option. You may learn about how to use Monte Carlo simulation 
to price exotic options in RMSC4001. Examples of exotic options may include asian option, binary 
option, barrier option, etc. 

Example 6.1. Forward start options: options that will start at some time in the future. 
Sometimes employee stock options (about corporate finance) can be viewed as forward start options. 
It is common to specify that the strike price will be set in the future so that the option is initially 
at the money or a certain percentage in the money or out of the money. 

Example 6.2. The payoffs from lookback options depend on the maximum or minimum asset 
price reached during the life of the option. Payoff: 

(i) Floating lookback call: max{0, St — miip St}. 

(ii) Floating lookback put: max{0, maxi St — St}- 

6.2 Risk-Neutral Valuation 

(a) Let II be the value of the portfolio consisting a long positing in A shares and a short position 
in one option (note that no specific option is assumed!). Therefore, at the end (time T) of the 
life of the option if there is an upward movement in the stock price, ILr = So'uA — v u . Similarly, 
if there is a downward movement in the stock price, Ily = SodA — v d . 

(b) To make the portfolio riskless, its value at maturity must be constant. That is, 

S 0 uA -v u = S 0 dA - v d . 

Rearranging the terms, we have 

* = Vu ~ Vd 

~ S 0 u ~S 0 d' 

(c) Since the portfolio is riskless, it must earn the risk-free rate r. Therefore, IIo = 1 ■ That 
is, 

SoA — v = (SquA — v u )e~ rT . 

Rearranging, we have 

v = SoA - (SquA - v u )e~ rT . (1) 
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6.3 Multi-period Binomial Trees 


6 APPENDIX 


(d) By direct substitution of A = — * — — into (1). We have, 

•Sq u — So d 


v = 


v u ~v d ( u{v u - v d ) 


u — d 


u — d 


- v n e 


—rT 


= e 


-rT 


rT v u -v d dv u — uv d 


u — d 


= e 


-rT 


e rT — d 
u — d 


V u + 1 ~ 


u — d 
e rT — d 


u — d 


v d 


e rT - d 

If we denote p* := , then we have 

u — d 


v = e rT \p*v u + (1 -p*)v d ]. 


If there is no arbitrage, then p* e (0, 1) (otherwise, buying a risk-free asset can always dominate the 
payoff of a stock or buying the stock can always dominate the payoff of a risk-free asset). Therefore, 
it is natural to interpret p* as the probability of the up movement. Under this probability, 

E *{S T ) = p* Sou + (1 -p*)S 0 d = p*S 0 (u -d) + Sod = S 0 e rT . 

Therefore, under this probability, the expected return on the stock equals the risk-free rate. This 
is an example of risk-neutral valuation and you will learn more in other RMSC courses. 

Remark 6.1. Note that no assumption on the payoff of the option is needed. That means, you 
can price any option in this way. 


6.3 Multi-period Binomial Trees 


With two periods, we have 


v u 

v d 

v 


e- rAt \p*v uu + (l-p*)v ud ] 

e- rAt \p*v du + (l-p*)v d d] 

e- rAt \p*v u + (l-p*)v d } 
e rAt - d 
u — d 


After substitution, we have 

v = e~' 2rAt [(p*) 2 v uu +p*( 1 -p*)v ud + p*( 1 -p*)v du + (1 -p*) 2 v dd \. 

Here p* is the risk- neutral probability. v uu is the payoff when u 2 S occurs. Other notations are 
defined similarly. For simplicity, I assumed u and d are the same at each node of the tree and so 
that p* is the same at each node. 

6.4 Examples 

Example 6.3 (Two-step Binomial Tree on Exotic Options). Consider a two-step binomial tree for 
the stock price S(t ) at t = 0, 1,2 year. The stock price starts at 5(0) = 100 and moves either up 
or down by 20% per year. Assume that the (continuously compounded) risk-free rate is 3% per 
annum. 
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6.4 Examples 


6 APPENDIX 


(a) Estimate the price of a European put option with strike K = 110 and expiry = 2 years. 

(b) Estimate the price of a European option with expiry at 2-year and the payoff function given 
by: 


Payoff = max{5(l) — 110,0} + max { 5(2) — 110,0} 
Solution, (a) First, calculate p*, 


P 


* 


e 0 03 - 0.8 
1.2 - 0.8 


0.57614. 


The payoffs of the put option when stock price is S uu , S ud and S dd are 0, 14 and 46 respectively. 
v = [0(p*) 2 + 14(2)(p*)(l - p*) + 46(1 - p*) 2 ]e -0 ' 03 ' 2 = 14.2226. 


(b) Note that v uu = 44, v ud = 10, v du = 0 and v dd = 0. 


v = [44(p*) 2 + 10(p*)(l — p*) -(- 0(p*)(l — p*) + 0(1 — p*) 2 ]e -0 ' 03 ' 2 = 16.0545. 
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